GL-TR-89-0279 

ENVIRONMENTAL  RESEARCH  PAPERS,  NO.  1041 


AD-A219  546 


A  Simple  Treatment  of  "Snowplow"  Models  of  Explosions 


Christopher  Sherman 


IONOSPHERIC  PHYSICS  DIVISION  PROJECT  231 0 

GEOPHYSICS  LABORATORY 


HANSCOM  AFB,  MA  01731-5000 


90  03  09  067 


'This  technical  report  has  been  reviewed  and  is  approved  for  publication' 


DAVID  N.  ANDERSON,  Acting  Branch  Chief 
Ionospheric  Modelling  &  Remote  Sensing  Branch 
Ionospheric  Physics  Division 


BOR  THE  COMMANDER 


RORERg^Z.  SKRIV.ANEK,  Director 
Ionospheric  Physics  Division 


This  report  has  been  reviewed  by  the  ESD  Public  Affairs  Office  (PA)  and  is 
releasable  to  the  National  Technical  Information  Service  (NT IS) . 


Qualified  requestors  may  obtain  additional  copies  from  the  Defense  Technical 
Information  Center.  All  others  should  apply  to  the  National  Technical 
Information  Service. 


If  your  address  has  changed,  or  if  you  wish  to  be  removed  from  the  mailing 
list,  or  if  the  addressee  is  no  longer  employed  by  your  organization,  please 
notify  GL/IMA,  Hanscom  AFB,  MA  01731.  This  will  assist  us  in  maintaining 
a  current  mailing  ist. 


Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  requires  that  it  be  returned. 


DD  FORM  1473,  JUN  86  Previous  editions  are  obsolete.  SECURITY  CLASSIFICATION  OF  this  PAGE 

UNCLASSIFIED 


Illustrations 


1.  Non-Dimensional  Kinetic  Energy  versus  Non-Dimensionalized  Expansion  4 

Radius  According  to  Present  and  Prev  us  Solutions 

2.  Non-Dimensionalized  Expansion  Radius  versus  Non-Dimensionalized  Time,  5 

According  to  Present  and  Previous  Solutions 

3.  Non-Dimensionalized  Expansion  Radius  versus  Non-Dimensionalized  Time,  6 

Including  Effects  of  Ambient  Pressure 


iv 


A  Simple  Treatment  of  “Snowplow” 
Models  of  Explosions 


1.  INTRODUCTION 

There  are,  aside  from  detailed  numerical  solutions,  several  approximate  ways  of  treating  the 
evolution  of  an  explosion  in  an  ambient  medium.  The  appropriate  approximation  for  a  given  case 
depends  on  the  nature  both  of  the  explosion  itself,  and  on  the  ambient  conditions  prior  to  the 
explosive  energy  release.  If  the  explosion  occurs  in  a  dense  medium,  the  whole  process  can  be 
described  by  continuum  mechanics,  and  a  shock  will  move  outward  into  the  initially  stagnant 
ambient  medium.  Here,  “dense"  means  that  the  radius  of  the  spherical  shi  ck  front  is  very  large 
compared  to  the  ambient  mean  free  path.  If  the  shock  is  also  very  strong,  the  pressure  outside  the 
shock  front  may  be  neglected  in  comparison  with  that  Inside;  this  reduces  the  number  of  initial 
parameters  entering  the  problem,  and  permits  solutions  of  a  special  type,  involving  ordinary,  rather 
than  partial  differential  equations.  Such  solutions  have  been  studied  in  some  detail.1-2  In  the 
opposite  limit,  namely  an  explosion  moving  into  a  vacuum,  there  can  be  no  exterior  shock.  This 
problem  has  been  studied  using  continuum  theory  for  an  initially  dense  interior  gas3  and  also  for 
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initially  tenuous  gases  having  a  Maxwellian  distribution,  using  a  free  molecular  flow  treatment.4 5 6 
These  two  approaches,  Interestingly,  give  quite  similar  results. 

It  is  for  cases  In  between  these  two  extremes,  when  the  ambient  density  is  neither  zero,  nor  dense 
enough  to  allow  the  use  of  continuum  theory  that  the  models  called  “snowplows"  apply.  Between  1964 
and  1969  several  articles  5>6'7  dealing  with  explosions  occurring  In  a  rarefied  atmosphere  appea.  d  In 
the  literature.  The  application  of  these  articles  was  primarily  to  high  altitude  chemical  releases  and 
all  three  utilized  a  more  or  less  detailed  hydrodynamical  treatment  of  the  expanding  gas  produced  by 
the  explosion.  This  involves,  in  general,  the  solution  of  a  set  of  non-linear,  partial  differential 
equations.  Here,  contrariwise,  the  internal  details  of  the  expansion  are  completely  ignored  and  only 
the  overall  conservation  of  momentum  is  considered.  It  is  the  purpose  of  this  note  to  illustrate  how  a 
simple  application  of  basic  principles  can,  on  occasion,  be  used  to  obtain  results  as  accurate  as 
extended,  much  more  complicated  calculations. 


2.  PRESSURE  FREE  EXPANSION 

Two  of  the  last  three  articles  cited  refer  to  a  “simple  snowplow"  model.  In  this  model,  it  is 
assumed  “that  the  sphere's  radius  expands  so  as  to  conserve  kinetic  energy,  T”.  Now,  the  expanding 
gas  that  results  from  the  explosion  is  considered  to  be  a  “spherical  piston",  which  as  it  expands, 
sweeps  up  ambient  gas  in  its  motion.  But  sweeping  up  gas,  that  is,  accretion,  is  clearly  an  inelastic 
process  and  it  cannot  be  correct  to  assume  that  kinetic  energy  is  conserved.  What  is,  to  a  first 
approximation,  conserved  is  outward  momentum;  and  the  conservation  of  momentum  may,  if  we 
assume  the  gas  to  expand  in  a  uniform  shell,  be  written  down  immediately  as 

M0Vo  =  MV  =  (M0  +  4/3  k  R3p)  V 

dR 

=  (M0  +  4/3  x  R3p)  ”  (1) 

Here,  t  is  the  time.  M0  is  the  mass  of  the  released  gas,  V0  the  initial  outward  velocity  of  the  shell,  M  and 
V  the  corresponding  quantities  at  later  times;  R  =  R(t)  the  expansion  radius:  and  p  the  density  of  the 
ambient  gas.  Introducing  a  non-dimensional  variable  <{>  given  by 

/Mn  \l/3 

R=(^J  *  <2> 

and  letting  T0  =  E0  be  the  Initial  kinetic  and  total  energy  respectively,  assumed  equal  at  t  =  O1,  we  have 


4.  Molmud,  P.  (1960)  Expansion  of  a  rarefied  gas  cloud  into  a  vacuum.  Phys.  Fluids  3:362. 

5.  Stuart,  G.W.  (1965)  Explosions  in  rarefied  atmospheres,  Phys.  Fluids  8:603. 

6.  Holway,  I.H.,  Jr.  (1969)  Similarity  model  of  an  explosion  in  a  rarefied  atmosphere,  Phys.  Fluids 

12:2506. 

7.  Klein,  M.M.  (1968)  Similarity  solution  for  cylindrical  gas  cloud  in  rarefied  atmosphere,  Phys. 

Fluids  11:964. 


2 


(3) 


T  =  1/2  MV2=  1/2  M 


M02V02 

M2 


=  1/2 


MQ2V02 

M 


T/T0 


Mo  Mn _ 1 

M  =  M0  +  4/3  JtR3p  “  1  +  «(>3/3 


(4) 


In  Figure  1,  we  plot  this  last,  superimposed  on  Figure  3  of  Holway,  and  see  immediately  that  this 
simple  procedure  gives  results  intermediate  between  those  of  Stuart  and  the  later  results  of  Holway. 

If,  instead  of  a  shell  expansion,  we  assume  a  similarity  solution,  v(r)  =  r/RV(R)  where  r  and  v  are 
the  radius  and  corresponding  velocity  at  points  interior  to  the  expansion  surface,  the  momentum  is 
given  by  3/ 4  MV;  and  the  kinetic  energy  by  3/5  MV2.  Eqs.  (  1)  and  (4)  remain  valid,  however,  as  does  the 
subsequent  integration.  It  is  assumed  in  this  case  that  the  interior  density  is  uniformly  distributed. 

A  natural  non-dimensionalization  of  the  time  leads  to 


t 


1_ 

V0 


'MqM/3  (  M05  M/6 

4npJ  S  =  [l28x2p2E03J 


s 


(5) 


The  last  form  is  identical  with  Holway’s  if  128  =  423(y-  l)3  or  y  =  5/3.  Equation  (1)  in  non- 
dimensional  variables  is 


(i  +  <>3/3)  ff-  =  1 


(6) 


or,  integrating,  with  <)>(0)  =  0, 

s  =  (J)  +  <{>4/12  (7) 

In  Figure  2,  Eq.  (7)  is  superimposed  on  Holway’s  Figure  1  and  again  gives  results  intermediate 
between  the  “simple  snowplow”  and  the  Stuart  solution  on  the  one  hand,  and  the  more  recent  solution 
of  Holway  on  the  other. 


3.  EXTERNAL  PRESSURE 

Although  all  three  of  the  referenced  articles  cany  their  solutions  up  to  and  well  beyond  values  of 
<S>  =  1,  none  of  them  includes  the  effects  of  an  extt  mal  pressure.  Insofar  as  these  models  are  applicable 
to  chemical  releases,  it  is  clear  however  that  for  such  radii  external  pressure  effects  must  play  a  role  in 
the  expansion;  and  In  fact,  that  for  expansion  much  beyond  =  1,  they  will  rapidly  come  to  dominate 
the  process.  This  will  occur  because  the  internal  momentum  is  fixed,  while  the  external  pressure 
exerts  a  force  which  is  proportional  to  the  expanding  surface  area. 

To  take  account  of  the  pressure  in  the  simplest  way,  we  replace  Eq.  (1)  by 

|r  [(Mo  +  4/3  x  R3p)  If  ]  =  -4x  R2  p„  (8) 

where  P0  is  the  ambient  pressure.  Now  ordinarily  p0  =  c2/3p.  (c  is  the  ambient  thermal  velocity.) 

Here,  however,  material  Is  not  reflected,  but  accretes;  hence  we  set  p0  =  c  2/6p  to  obtain 
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Figure  1.  Non-Dimensional  Kinetic  Energy  Versus  Non-Dimensionalized  Expansion  Radius  According  to  Present  and 
Previous  Solutions.  Here  l  is  related  to  an  (assumed)  form  of  the  density  distribution  (uniform,  shell,  etc.)  and  y  is  the  ratio 
of  specific  heats. 


CD  SIMPLE  SNOWPLOW  (STUART) 
©  STUARTS  SOLUTION 


■Dimensionalized  Expansion  Radius  Versus  Non-Dimensionalized  Time,  According  to  Present  and  Previous  Solutions. 


is  the  initial  expansion  Mach  number. 


(9) 


+  4/3  n  R3p)  4rcR2p(^)  2=  -4tcR2P  c2/6 


Non-dimensionalizing,  as  before 


Letting  a  =  sonic  velocity  of  the  ambient  gas,  a2=  2  so  that 

(l  +  <t>3/3)  ~|+  4>2  (  ff-)  2=  -<1>2/2yM2  (11 

where  M,  the  initial  (“hybrid")  Mach  number  will  not  be  confused  with  total  mass.  A  first  integral  of 
Eq.  (11)  is 


(  d£\  2  _ 

Vds)  ~  (i  + 


C  1 

(1  +  <)>3/3)2  "  2yM2 


When  t  =  0.  <t»  =  0.  and  |o=l  sothatC=l  +  ~2  30(1 

fd4^2  i  l~  2y  M  2  +  1  1 

Ids J  ~  2yM2  [(1  +  4> 3/3)2  _1J 


(13) 


An  absolute  upper  limit  on  (J)  is  now  given  by  ( 1  +  <j>3/3)2  <  2y  M  2  +  1  or  (2  +  4>3/3)  <  67  M  2 .  If  7  =  5/3 
and  M  =  10  (typical  values),  tlien  <J>  <  2.5;  a  value  smaller  than  some  of  those  plotted  by  Stuart  and 
Holway.  Because  of  the  <\>6  term,  this  result  is  rather  insensitive  to  the  value  of  M. 

Taking  the  root  of  Eq.  (13)  and  separating, 

V27M(1  +  d>3/3)  d<}>  _  (14 

[2y M 2  +  1  -  (1  +  4>3/3)2] */2 


The  left  hand  side  of  Eq.  (14)  leads  to  an  elliptic  integral  and  cannot  be  integrated  In  elementary  terms. 
It  is  most  conveniently  handled  by  a  numerical  integration,  and  in  Figure  3,  taking  7  =  5/3,  we  show 
the  results  of  the  integration  for  several  values  of  the  initial  Mach  number,  M.  Since  M  =  °o  is 
equivalent  to  the  pressure  free  case,  the  inclusion  of  the  external  pressure  brings  the  values  of  4>  closer 
to  those  obtained  by  Holway;  although  it  must  be  noted  that  the  validity  of  snowplow  solutions  for 
values  of  4>  >  1  is  questionable. 

In  spite  of  the  extreme  simplicity  of  these  considerations,  the  results  obtained,  both  with  and 
without  external  pressure,  compare  favorably  with  those  of  the  earlier,  more  detailed  calculations. 
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